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00 ' Abstract 

, The ^^F production during the first dozen thermal pulses of AGB stars with masses 

^ ■ M and metallicities Z (M = SMq, Z = 0.02), {M = GM©, Z = 0.02) and {M = 3Mq, 

I Z = 0.001) is investigated on grounds of detailed stellar models and of revised rates for 

■ -'^^N(a,7) ^^F and ^®0(q;, 7)^^Ne. These calculations confirm an early expectation that 
I -"^^F is produced in AGB thermal pulses. They also enlarge substantially these previous 
' results by showing that the variations of the level of ^^F production during the evolution 

is very sensitive to the maximum temperature reached at the base of the pulse. These 

Q ' variations are analyzed in detail, and are shown to result from a subtle balance between 

I different nuclear effects (mainly ^^F production or destruction in a pulse, and ^^N synthesis 

^2 ' during the interpulse), possibly superimposed on dilution effects in more or less extended 

I pulse convective tongues. 

■ Our calculations, as most others, do not predict the third dredge- up self-consistently. 
k>( , When parametrized, it appears that our models of intermediate-mass AGB stars are able 

}_j ■ to account only for the lowest ^^F overabundances observed in solar-metallicity MS, S 

I and C stars. That conclusion is expected to hold true for low-mass stars when fluorine 



production results from secondary ^■^C only. Massive AGB stars, on the other hand, 
are not expected to build up large surface F abundances. Therefore, the large fluorine 
overabundance reported for the super Li-rich star WZ Gas (where hot bottom burning 
is supposed to be operating) remains unexplained so far. Our results for the (3Mq, 
Z = 0.001) star indicate that fluorine surface overabundances can also be expected in 
low-metallicity stars provided that third dredge-ups occur after the early cool pulses. 
The relative increase in the surface ^^F/^^C ratio is, however, lower in the low-metallicity 
than in the solar-metallicity star. No observations are reported yet for these stars, and 
are urgently called for. 
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1 Introduction 



Fluorine overabundances in giant stars of spectral types MS, S and C have been reported 
recently (Jorissen, Smith & Lambert 1992; Paper I). These observations reveal a progressive F 
enrichment of the envelope as the C/0 ratio increases, reaching [F/0] ~ 1.5 in some C stars. 
An independent indication of large F overabundances in carbon stars comes from the recent 
detection of AlF in the inner envelope of the dust-enshrouded carbon star IRC+ 10216 (Ziurys, 
Apponi & Phillips 1994). These observations shed light on the much debated question of the 
nucleosynthetic origin of ^^F. They clearly point towards thermal pulses as a likely site for 
its production. 

The first quantitative proof that ^^F can indeed be produced in thermally pulsing AGB 
stars has been provided by Forestini et al. (1992; Paper II). However, this early investigation 
just considered the first 4 thermal pulses of a single (3 Mq, Z = 0.03) star. As a conse- 
quence, it left unanswered several important questions, like the dependence of the level of 
I9p production on the stellar mass, metallicity, or pulse number. 

The present paper aims at exploring some of these questions. It analyzes a dozen pulses in 
three stars of different masses and metallicities (M = 3Mq, Z = 0.02; M = OMq, Z = 0.02; 
M = 3M0, Z = 0.001). Use is made of an improved code which follows more accurately the 
structural evolution of the AGB stars and the concomitant nucleosynthesis. It also updates 
some nuclear reaction data of importance in the study of the ^^F production. In particular, 
a more accurate rate (de Olivera et al. 1995) for the key reaction ^^N(a,7)^^F is adopted, 
as well as a recent re-evaluation (Giesen et al. 1994) of the ^^0(a,7)^^Ne rate. In addition, 
the surface fluorine enhancement due to the first and second dredge-ups is also investigated. 

The nuclear transformations leading to the production of ^^F in thermal pulses are re- 
viewed in Sect. |2|. The evolution code and its key ingredients are briefly described in Sect. ^, 
while Sect. |^ presents our calculated ^^F abundances. Our expectations are confronted with 
the observational data in Sect. ||, and conclusions are drawn in Sect. ^ 

2 The nuclear transformations governing the ^^F production 
in thermal pulses 

The reaction chain producing ^^F in thermal pulses has been identified in Paper II, and is 
displayed in Fig. ||. Its main characteristics may be summarized as follows: 

(1) the main seeds for ^^F production are ^^C and ^^N left behind by the H-burning shell. 
The resulting relatively small amount of ^^C is the main factor that limits the amount of ^^F 
possibly produced in a thermal pulse. The efficiency of the ^^F synthesis in a pulse can be 
expressed as (yo^t(^^N)+Y*^t(^^F))/Yin(^^C). In this ratio, ^^("'^^C) represents the abundance 
by number of ^^C injected in the pulse, and yoy^(^^N) and ^^^^^(^^F) are the abundances by 
number of ^^N and ^^F at the end of the pulse, assuming zero initial ^^F and ^^N abundances. 

(2) the main loss of efficiency in the ^^F synthesis comes from the (n,7) reactions, mainly 
on °"Fe and heavier nuclei. They limit the occurrence of ^4N(n,p)i4C and 26A1 (n , p) 26Mg 
reactions, which are the main producers of the protons required for manufacturing ^^N, the 
^^F progenitor]^. 

^The role of ^®A1 in enhancing the ^^F production in thermal pulses has been identified in Paper II. The 
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Figure 1: The chain of reactions producing ^^F in thermal pulse conditions. Bold italicized 
species are ashes from the H-burning shell. The most important (n,7) capture competing 
with the indicated (n,p) reactions is ^^Fe (n , 7) 

(3) almost all the protons liberated as indicated above can lead to the production of ^^F if 
indeed the ^^O (p , a) ^^N (a, 7)^^F chain develops. A necessary condition for this to occur is 
the availability of ^^O. Some ^^O might be inherited from previous pulses. As no interpulse 
^^0 burning is predicted by our models, this requires mainly that the duration Atpuise of these 
pulses be shorter than the destruction time scale of ^^O against a-capture. If this is not the 
case, no ^^O from former pulses can survive at the beginning of a pulse. In such conditions, 
and as ""^^N (a , 7) ""^^F {(3~^) ^^O is slower than -"^^C (a , n) ^^O for typical pulse temperatures, 
one would expect the neutrons to be liberated and captured by the heavy nuclides before any 
^^0 is made available. In this case, the reaction chain leading the synthesis of ^^F displayed 
in Fig. 1 would not operate. 

In fact, the situation as it comes out of the detailed calculations to be reported in Sect. ^ 
is quite different. In reality, the nuclear sequence of Fig. |l] can still operate in a pulse even if 
that pulse does not inherit any ^^O from previous pulses. We refer to Sect. |^ for details. 

(4) ^^F has the rather large Maxwellian-averaged neutron-capture cross section of 5.7 mb 
at 30 keV (Beer, Voss & Winters 1992). It has been checked that (i) ^^F(n,7)2°F has no 
influence on the neutron budget. This results from the fact that ^^N or ^^Fe are much more 
important neutron captors than ^^F; and (ii) ^^F (n , 7) ^'^F could become a significant ^^F 
destruction channel only for much larger neutron fluences than the ones predicted by the 
models reported in this paper. In fact, if the larger neutron production would result from a 
primary ^^C reservoir, an increased ^^F production would also be obtained from the chain of 
transformations of Fig. 1, thus limiting the impact of ^^F (n , 7) ^'^F . The models reported 
in this paper do not consider such a possibility, which deserves further detailed studies. 

(5) as shown in Fig. 1, ^^F (a ,p) ^^Ne could destroy some of the produced ^^F if the corre- 
sponding destruction time scale, TaC^F), gets shorter than Atpuise- This happens when the 
temperature at the base of the pulse exceeds typically ~ 300 10^ K. The ^^F production is 
therefore expected to be less efficient in the late AGB evolutionary stages, when the pulse 
temperatures are high. 

H-burning shell of AGB staxs is known to produce large amounts of ^"^Al (Forestini et al. 1991) 
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Figure 2: The 6th and 7th pulses of the 3 Mq, Z = 0.02 star. The scales are the same 
in both panels. Hatched regions correspond to convective zones. Each computed model 
is represented by a vertical line in the convective zones (horizontal lines have no special 
meaning) . The short-dashed lines identify the location of maximum energy production in the 
H-burning (top) and He-burning (bottom) layers. Regions a and h as defined in the text are 
also indicated 



Finally, it has to be noted that the F synthesis in convective pulses involves species that 
may have nuclear lifetimes shorter than, or similar to, the convective mixing time scale. This 
situation is encountered for neutrons and protons. A proper treatment of the ^^F produc- 
tion may therefore require the use of an algorithm coupling nucleosynthesis and convective 
diffusion. This question will be addressed further in Sect. ^. 



3 Key ingredients of the stellar evolution code 

A new implicit stellar evolution code that differs in several ways from the one used in Paper H 
has been designed in order to follow accurately the evolution of the structurally very complex 
AGB stars (see Mowlavi 1995a for a full description of the code). Its new features significantly 
improve the numerical accuracy of the models, and avoid spurious numerical effects exhibited 
by the models of Paper H. 

In particular, the algorithm defining the discrete mass zones is based on the structural 
equations, rather than on the relative variations of the dependent variables. This provides 
a better handle on the accuracy achieved on the structure]^. Convective zones growing into 
regions of variable chemical composition have been modeled with great care. A special 
algorithm moreover ensures that composition discontinuities remain sharp as the star evolves, 
thus avoiding "numerical chemical diffusion" as encountered in the models of Paper II. 

Each model is computed by iterating the nucleosynthesis -|- mixing and the structure 
calculations. Starting with an initial model at time t, taken identical to the model at time 

^This structural accuracy is different from the relative precision required for the convergence of the Newton- 
Raphson solution of the structural equations. The latter easily reaches 10~^ in all our models. 
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t — At, where At is the time step between the two models, a first iteration is performed by 
calculating the nucleosynthesis + mixing, followed by the calculation of the structure through 
a Newton-Raphson method. A second iteration performs the nucleosynthesis + mixing on 
the last calculated structure, followed by a new calculation of the structure. The procedure 
goes on until the structure and the chemical distribution have converged. Two to three 
such iterations are usually sufficient. From the core helium-burning phase on, however, the 
number of iterations is limited to two, still providing a satisfactory accuracy within reasonable 
computer times. Rezoning is applied after each structure calculation in order to reach the 
required accuracy, and to predict the location of the convective boundaries to better than 10~^ 
M.Q. Instantaneous mixing and concomitant homogeneous composition are assumed in the 
convective zones, except for nuclides with nuclear lifetimes Tnuc shorter than the convective 
mixing time scale Tmix- More specifically, the nucleosynthesis is calculated in a one-zone 
approximation making use of an effective reaction rate obtained by averaging the reaction 
rate over the convective zone (see Prantzos, Arnould & Arcoragi 1987 for more details). When 
■^nuc < '^mix for a given nuclide, its abundance must be computed separately in each layer. 
In practice, it is found that these nuclides (mainly neutrons, protons and sometimes ^^O) 
usually also satisfy the condition of local equilibrium between production and destruction, 
from which their abundance can easily be computed in each layer. This procedure is only 
applied to nuclides with X < 10^^^ in the convective zone (where X is the mass fraction), 
like neutrons and often protons. Species more abundant than X — 10 or with t^uc ^ '^mix 
require an algorithm coupling nucleosynthesis and convective diffusion. Such an algorithm 
has been implemented in our code, but is very time-consuming. It has been applied to the 
calculation of one pulse, as described in Sect. 4T, in order to evaluate the accuracy of the 
results derived from the simpler procedure. 

Nuclear reaction rates are taken from Caughlan & Fowler (1988; CF88), except for the 
reactions listed in Appendix. New reaction rates of importance for the ^^F nucleosynthesis 
include ^^N{a,-fy^F (de Oliveira et al. 1995) and ^^0{a,'y)'^'^'Ne (Giesen et al. 1994). The 
new ^^N(a,7)^^F rate is slower than the CF88 one by a factor 1.5 to 50 in the temperature 
range 2 < Tg < 3 of interest here (Tg is the temperature in units of 10* K). The impact of 
this change on the ^^F nucleosynthesis has been investigated by de Oliveira et al. (1995). 



4 Model predictions 

Three intermediate-mass stars are computed with (mass, metallicity) pairs of (3, 0.02), (3, 
0.001) and (6, 0.02), referred to as M3Z02, M3Z001 and M6Z02, respectively. 

Our models are computed from the pre-main sequence phase, when the star is on the 
Hayashi track with central temperatures less than 10® K. The structural and chemical evo- 
lution of the three stars are followed all the way up to the 13th, 15th and 11th pulse of the 
AGB phase for the M3Z02, M3Z001 and M6Z02 stars, respectively. A detailed description 
of the structural and chemical evolution of these stars can be found in Mowlavi (1995a). In 
the remaining part of this paper, we concentrate only on the fluorine production in the AGB 
phase. 
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Figure 3: Abundance profiles before (left panel), and after (right panel) the 7th pulse of the 
3 Mq, Z=0.02 star. Regions a and b defined in the text are readily identified on the left panel. 
In the left panel, the rapid drop of the ^^F abundance to the right of the plateau corresponds 
to the maximum outward extension of pulse 6, and marks the transition between regions a 
and b. The dots on the ^^F profiles correspond to the mesh distribution 



4.1 The 3 Mq, Z = 0.02 case 

The features governing the ^^F synthesis in relatively cool thermal pulses can be illustrated 
by considering the 6th and 7th pulses of the M3Z02 star. The location of the convective 
boundaries of these pulses is displayed in Fig. |2| as a function of time. The pulse duration, 
Atpuisei is 182 and 159 y for the 6th and 7th pulse, respectively, and the maximum base 
temperature Tpuise,b = 239 and 245 10^ K, respectively. Defining the maximum inward and 
outward extensions of pulse n as M^"jj and Mp^^ , respectively, the overlap 

between pulses 6 and 7 amounts to 0.753. Thus 75.3% of the mass of pulse 7 (defined as 
region b on Fig. ^, and extending from Mr = 0.572 to 0.591 M©) is inherited from pulse 6, 
the remaining fraction (defined as region a on Fig. ^, and extending from Mr = 0.591 to 
0.597 Mq) containing the ashes left behind by the hydrogen-burning shell. These two regions 
are readily identified on the left panel of Fig. ^, which displays the abundance profiles of ^^C, 
i^N, i^N, 1^0, ^^F and ^^Al before pulse 7. The drop in ^^N and ^^O observed at the bottom 
of region b translates the fact that a-capture reactions continue to operate on these nuclides 
in the hottest layers below pulse 6. A concomitant small hump in the fluorine abundance is 
also observed in these layers. On the other hand, the non-homogeneous proflle exhibited by 
^^O in the upper part of region b corresponds to the ^^O abundance frozen, in a given layer, 
at its value in the 6th pulse at the time the considered layer became radiative. 
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Figure 4: Evolution of the ^^N, ""^^N, ^^O and ""^^F mass fractions during the 7th pulse of the 
3M0, Z=0.02 star. The evolution of the temperature at the base of the pulse is also shown 
(solid line), to be read on the right scale (in units of 10® K). The quantities are plotted every 
4 to 6 models, resulting in a poor time-resolution for some species 



The evolution of the abundances of several nuclides of interest during the 7th pulse is 
shown in Fig. ^ The ingestion of ^^N and the concomitant production of ^^O are clearly 
visible during the growing phase of the pulse. The ^^F abundance, on the other hand, first 
decreases slightly as the pulse grows into fluorine-depleted regions. It then increases as soon 
as ^^N (a , 7) ^^F starts operating. A net increase in the fluorine abundance at the end of the 
pulse, as compared to its value before the pulse, is clearly apparent on Figs. ^ and ^, even 
though the conversion of ^^N into ^^F has not been complete. It is noteworthy that such a net 
^^F production has been possible even though the amount of ^'^O inherited from the previous 
pulses is about ten times lower than the ^^F abundance level. In fact, ^^C and ^^N are not 
ingested instantaneously at the start of the pulse, but are instead mixed in progressively. The 
ingestion time scale is set by the rate of outward progression of the pulse convective tongue. 
This time scale turns out to be longer than the one for ^^O production. In other words, 
neutrons are released by ^^C (a , n) in presence of ^^O produced earlier in the pulse by 
i^N {a , 7) (/3+)i80 (see Fig. |). 

The intershell mass fractions Xout of ^^N and ^^F at the end of each pulse are shown 
in Fig. |5| as a function of the pulse number. The fraction of ^^N converted into ^^F is low 
in the first pulses ^. It is, however, an increasing function of the pulse number, i.e. of the 
maximum temperature reached at the base of the pulse. The dotted line in Fig. ^ represents 

^The conversion of ^^N into ^^F was faster with the CF88 rate for ^^N (a , 7) ^^F . We refer to de Ohveira 
et al. (1995) for a discussion of the impact of the new rate on ^^F production in thermal pulses 
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Figure 5: Various quantities in the SMq, Z=0.02 star as a function of pulse number. Lower 
panel: Maximum temperature at the base of the thermal pulse; Middle panel: Intershell ^^N 
(filled circles), ""^^F (filled squares) and ^^N+^'^^F (dotted line) mass fractions at the end of 
the pulse, as well as ^^C mass fraction averaged over region a as defined in Fig. |2| (open 
triangles); Upper panel: Efficiency factor g of fluorine production in the pulse (filled squares) 
and dilution factor of the hydrogen-burning ashes into the pulse (open circles) 



the maximum ^^F mass fraction that would result from the complete conversion of ^^N into 
^^F. That limiting value is roughly constant with pulse number. In other words, the efficiency 
of ^^N production through the above-mentioned reactions is not very sensitive to the pulse 
temperature, at least when Tpuise.b < 260 10^ K. The fraction of ^^N that can be converted 
into ^^F in these pulses is however sensitive to the pulse temperature (Fig. ^). 

As already indicated in Sect. |2|, the amount of ^^N that can be produced is limited by 
the available neutrons, i.e. by the available ^^C supply. The maximum efficiency would 
correspond to a situation where each neutron liberated by ^'^C {a , n) ^^O produces a proton 
which is then captured by ^^O (p , a) ^^N . The actual efficiency can be derived from baryon 
number conservation, as expressed by the relation 

_ ^|Xo,t(^^N) + Xo,t(^^F) - iXin(^5N) - XiX^F) 

where the subscript in refers to the intershell mass fractions averaged over region a + b before 
the considered pulse, and the subscript out to the abundance at the end of the pulse. The 
average ^^C mass fraction injected in the pulse, Xin(^^C), results from the dilution into the 
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pulse of ^'^C left behind by the hydrogen-burning shell [Xa(^^C) in the middle panel of Fig. ^. 
The dilution factor 

ra = Ma/{Ma+Mh), (2) 

where Ma and Mf, are the masses of regions a and b, respectively (see Fig. amounts to 
about 25% for the 7th pulse. The actual efficiency turns out to be around 20%, as seen in 
the upper panel of Fig. ^. 

Another way of evaluating the efficiency is based on nuclear flux considerations. Leakage 
on the path from ^^C to ^^N first occurs from neutrons liberated by ^'^C (a , n) ^^O and 
involved in (n,7) rather than in (n,p) reactions. The main neutron poison in that respect is 
^^Fe. The efficiency gi with which protons are liberated from a ^^C seed is 

E <XA,npi1n(^Z) 
n— capture 

CTl4.nplin('^N) + a26,r.pYin(''Al) 
C7i4,np>in(l^N) + dae.npYin (26A1) + CT56,„Yi„ (S^Fe) ' 

where Y{^Z) = X(^Z)/A and a stands for the Maxwellian-averaged neutron-capture cross 
section. The sum E* i"uns over all n-capture reactions except ^^C (n , 7) ^'^C , as this reaction 
replenishes ^^C. Since ^"^N is only present in the hydrogen-burning ashes (region a; see Fig. 
its abundance after dilution in the pulse is Xin(i'^N) = raXa(^''N), with Xai^^N) = IQ-^, Xa 
designating the mass fraction averaged over region a. A similar relation holds for ^^Al in the 
M3Z02 star, with Xa(^'^Al) ~ 510"^. On the contrary, Xin(^'^Fe) is independent of ra, as 
^^Fe keeps its primordial abundance throughout the star. Therefore, as ^^Fe constitutes the 
main neutron poison competing with ^^N(n, p)^'^C and ^^Al(n, p)^^Mg, its importance will be 
comparatively larger, and thus gi smaller, for smaller dilution factors r^. The efficiency of 
fluorine production in the pulse is thus dependent on the pulse overlap factor r = 1 — ra, gi 
increasing with decreasing r. This is clearly seen on the upper panel of Fig. ^, which shows 
that the efficiency correlates well with the dilution factor. On the contrary, the efficiency is 
seen to be almost independent of the temperature, at least over the range 1.8 < Tg < 2.5. 
This is consistent with the fact that the Maxwellian-averaged rates for the neutron-capture 
reactions of interest are almost independent of temperature. 

A second leakage in the fluorine production comes from protons captured by nuclides 
other than ^^O. The efficiency 52 with which protons liberated by (n,p) reactions contribute 
to ""^^O (p , a) ^^N may be written 

E ^A,p^1n(^Z)' 
p— capture 

where the symbols are as in Eq. (^). The analysis of the reaction fluxes indicates that 
g2 ~ 0.90 in the pulses of the M3Z02 star. 

Finally, the efficiency of fluorine production can be decreased further by ^^F (n , 7) ^'^F 
and/or ^^F (a , p) ^^Ne . The n-capture reaction does not play an important role given the 
low neutron fluences encountered in our calculations. The a-capture reaction, on the other 
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hand, is not (yet) operating in the (stih too cold ) pulses encountered in the M3Z02 star. If 
(1 — (73) is the fraction of ^^F destroyed, the total efficiency g of fluorine production is given 
by 5 = 91 * 52 *53- 

The maximum fluorine abundance that can be expected in the intershell region can now 
easily be estimated from the above considerations. Assuming that ^^N is fully converted into 
^^F during the pulses (as is the case after a few pulses), the ^^F mass fraction at the end of 
pulse n is equal to 

X^X'^) = (1 - ra)xi°-^)(i9F) + r,i^gXH(i='C), (5) 

where the first term on the right-hand side represents the fluorine abundance inherited from 
the previous pulses, and the second term the amount of fluorine produced from the ^^C 
injected in the pulse. As fluorine builds up pulse after pulse, its abundance will reach an 
asymptotic value when X^^^ = X^^^ (assuming a constant overlap factor). Equation (^ 
yields, in the asymptotic regime, 

^out(^'F) = ^gX,(i3c) (6) 

This value can easily be computed from the equilibrium ^^C mass fraction XaQ^^C) resulting 
from the CNO cycles, and from the neutron balance expressed by Eq. (^). With Xai^^C) ~ 
10~^ and g ~ 0.20, this leads to Xout(^^F) — 310~^, which is precisely the asymptotic value 
obtained in the M3Z02 star (Fig. ||^. 

The results presented above are obtained assuming instantaneous mixing in the convective 
zone (except for neutrons). This assumption is, however, not strictly valid for protons and 
in some cases ^^O, whose nuclear lifetimes may become shorter than the convective mixing 
time scale. The use of a very computer time-consuming algorithm coupling diffusion and 
nucleosynthesis (Sect. ^ does not lead, however, to significantly different ^^F abundances in 
a test run on the 7th pulse. 

4.2 The 3 M©, Z = 0.001 case 

The evolution of the fluorine abundance in the SM©, Z=0.001 star is shown in Fig. |^ as a 
function of the pulse number. Several important differences exist as compared to the M3Z02 
case: 

i) The maximum fluorine abundance, corresponding to the complete transformation of ^^N 
into ^^F, is reached much faster in the M3Z001 star (after 3 pulses) than in the M3Z02 star 
(at least 10 pulses are needed). This is due to the higher temperatures reached in the pulses 
of the low-metallicity star. 

ii) The efficiency of ^^F synthesis (upper panel of Fig. ^) during the first pulses is much larger 
than in the M3Z02 star {g ~ 0.50 as compared to 0.20). This is due to the combined effect 
of smaller overlap factors between successive pulses (see the discussion of Sect. |4.1| ) and of 
larger ^^Al/^^Fe abundance ratios left behind by the hydrogen-burning shell and ingested by 
the pulse. Arnould & Mowlavi (1993, see also Arnould et al. 1995) showed that the ^^Al 
abundance resulting from hydrogen burning sensitively depends upon the temperature, being 
maximum for temperatures around Tg ~ 70. The temperature Thbs in the hydrogen-burning 
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Figure 6: Same as Fig. ^, but for the 3Mq, Z=0.001 star. In the middle panel, Xcno(^^C) 
corresponds to the ^^C mass fraction left behind by the hydrogen-burning shell. The difference 
between Xcno(^'^C) and Xa(^^C) is due to the interpulse ^^C burning 

shell is indeed more favorable to ^^Al production in the M3Z001 star (Thbs = 73 lO*' K 
before the first pulse) than in the M3Z02 star (Thbs = 45 10^ K before the first pulse). As a 
consequence, the ability of ^^N and ^^Al to compete with the neutron poisons is much higher 
during the first pulses of the M3Z001 star. In the second pulse of the M3Z001 star, about 
10 and 50% of the neutrons liberated by ^'^C (a , n) ^^O are captured by ^''N (n , p) ^^C and 
^^Al (n , p) ^^Mg , respectively. 

The efficiency, however, rapidly decreases, and drops to zero after the 5th pulse. This 
results from two effects. First, the intershell ■^^Al abundance decreases as the temperature of 
the hydrogen-burning shell exceeds 70 10^ K. Secondly, the temperature in the He-rich layers 
is high enough for ^^C (a , n) ^^O to destroy ^'^C in the deep layers of region a during the 
interpulse period. This is clearly seen on Fig. |^, which displays the abundance profiles before 
the 5th pulse of the M3Z001 star. As a consequence, the surviving ^^C is injected in the pulse 
at a time where the base temperature is close to its maximum. At these high temperatures, 
the destruction time scale of ^^N against a-capture becomes smaller than its injection rate in 
the pulse, so that its abundance decreases. Both effects drastically reduce the ability of ^^N 
and ^^Al to compete with ^^Fe for capturing neutrons, resulting in a drop of the efficiency of 
fluorine production within the thermal pulses. 

iii) The interpulse ^'^C burning leads to concomitant ^^N production in radiative layers, as 
seen on Fig. ^. When injected in the pulse, this ^^N will directly produce ^^F through 
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Figure 7: Same as Fig. 0, but for the 5th pulse of the SM©, Z=0.001 star 



^^N (a , 7) ^^F . Figure ^ shows how the different contributions to the after-pulse fluorine 
abundance vary with pulse number. In the early pulses, ^^N is produced from ^^C within the 
pulse, with an efficiency decreasing with the pulse number. The after-pulse fluorine abun- 
dance, resulting from the conversion of ^^N, reaches a maximum at pulse 3 and then decreases. 
From pulse 5 on, a substantial amount of ^^N is produced in region a during the interpulse 
period [Xin,a(^^N)], and represents, after injection in the pulse [Xin(^^N)], the main contri- 
bution to the after-pulse fluorine abundance. It is noteworthy that the F plateau reached 
with this radiative ^^N contribution is lower than the one that would have resulted had ^^N 
been synthesized in the pulse. This translates the fact that the efficiency grad 

(- 0.09) of 

production in radiative layers is actually lower than in convective (cold) pulses. Indeed, in 
a radiative layer, the ^^O necessary to produce ^^N through ^^O (p , a) ^^N becomes avail- 
able through -^^N (a , 7) ^^F {P~^)^^0, only after a substantial amount of ^^C has already been 
burned. This results from the longer ^'^N {a ,^)^^F {P~^y^O time scale as compared to the 
i3C(a,n)i60 one. 

iv) Finally, fluorine destruction through ^^F (a , p) ^^Ne begins to operate when the temper- 
ature at the base of the pulse exceeds 300 10® K. This happens from the 13th pulse on in the 
M3Z001 star, and is clearly visible in Fig. ^, which shows the time evolution of the abun- 
dances in the 13th pulse. The efficiency of fluorine production decreases to negative values as 
a result of the fluorine destruction in the hot pulses of the low-metallicity star (Fig. ^, upper 
panel). The concomitant decrease of the after-pulse ^^F abundance is also observed in Fig. ^ 
(middle panel). 

4.3 The 6 Mq, Z = 0.02 case 

The situation for the M6Z02 star is qualitatively similar to that described for the M3Z001 
star. The evolution of the after-pulse ^^F is simply delayed in the M6Z02 star as compared 
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Figure 8: Mass fractions X^^, averaged over the total mass covered by the pulse, of ^^C (open 
triangles) and ^^N (open circles, long-dashed line) before each pulse of the SMq, Z=0.001 
star. The contributions from regions a and 6 (see Fig. |2|), to the ^^N injected in the pulse are 
represented by open circles on doted lines, as labeled in the figure. The fluorine mass fraction 
Xout after each pulse, averaged over the mass of the pulse, is given by the filled squares 



to the M3Z001 star, as a result of the lower pulse temperatures (namely Tq = 260 and 290 
in the 8th pulse of the M6Z02 and M3Z001 stars, respectively). 

A specific feature of the M6Z02 star concerns however the temperatures at the base of 
the convective envelope, Tenv.b, which are much higher in the M6Z02 than in the M3Z02 and 
M3Z001 stars. At the first pulse, Tenv.b amounts to 2.5, 6 and 50 10^ K in the M3Z02, M3Z001 
and M6Z02 stars, respectively, and to 4.6, 25 and 79 10^ K at the 10th pulse, respectively. 
As a consequence, hydrogen burning occurs in the envelope of the M6Z02 star. This feature 
is known as hot-bottom burning (HBB). 

The effect of hydrogen burning on fluorine depends on the temperature, as shown by 
Arnould Sz Mowlavi (1993). For Tq < 20, fluorine is produced by ^^O (p , 7) ^^F , by a factor of 
up to 3 with respect to its solar-system abundance. At higher temperatures, it is essentially 
destroyed by ^^F (p , a) ^^O . As a consequence of HBB, the ^^F abundance in the envelope of 
the M6Z02 star is thus expected to first increase, and then to decrease as the temperatures get 
higher. Figure 11 shows that the rna-ximurn envelope ^^F cibundajnce is rea^clied when 7ej^v,b — 
7010*5 K[|, though the overabundance remains smah {[^^F]=log{Xenv{^^F)/XQ(^^F)) = 0.07). 
More important is the fluorine destruction in the envelope as the temperature exceeds Tq = 70, 
which happens from the 4th pulse on in the M6Z02 star. Some fraction of the fluorine brought 
in the envelope by the third dredge-up (which we did not obtain in a self-consistent way in 
our models; see Sect. |5.2|) would thus be destroyed as a result of the HBB. 



^The "efFective" temperature characterizing the nucleosynthesis in a convective zone is lower than the 
temperature at the base of that zone, because temperature-dependent reaction rates must be averaged over 
the whole convective region (see Sect. ^) 
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Figure 9: Same as Fig. ^, but for the 13th pulse of the SMq, Z=0.001 star 



5 Comparison with observations 



Observations of fluorine abundances at the surface of giant stars are available only for stars 



with metallicities close to solar (Paper I), as displayed on Figs. 13 and 14. We thus restrict 
the comparison with observations to the M3Z02 and M6Z02 stars. Predictions for the low- 



metallicity star are nevertheless briefly presented in Sect. 5.3 



5.1 First and second dredge-ups 

Before considering the fluorine production during the thermally-pulsing AGB phase, we ad- 
dress the question of surface ^^F alterations resulting from the first and second dredge-ups. 
They occur shortly after the end of the central hydrogen- and helium-burning phases, re- 
spectively, when the envelope engulfs the deep layers processed by hydrogen burning (the 
operation of the second dredge-up being restricted to stars with M ^ 4M0). 

The possibility that the surface fluorine abundance may be altered by the first or second 
dredge-ups is suggested by the observation of otherwise normal K and M giants exhibiting 
i9py'i6Q ];-atios significantly larger than solar (Fig. |l3| ). Since these stars do not bear the 
typical signatures of AGB nucleosynthesis (e.g. overabundances of carbon and elements 
heavier than iron), fiuorine production in the thermal pulses on the AGB cannot be invoked 
to account for their larger-than-solar ^^F abundances. 

As mentioned in Sect. |4.3| , fluorine is expected to be produced as well in hydrogen-burning 
layers when Tg < 20 (Arnould & Mowlavi 1993, Arnould et al. 1995), and destroyed at higher 
temperatures. This gives rise to a peak in the ^^F profile, as shown in Fig. |l^. The first 
dredge-up leads to a moderate increase (by about 0.05 dex) of the surface ^^F abundance in 
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Figure 10: Same as Fig. |, but for the 6M0, Z=0.02 star 



the M3Z02 star, along with a decrease of the ratio (Fig. |T^). No second dredge-up 

is experienced by this star. In the M6Z02 star instead, the temperatures in the hydrogen- 
burning layers are higher, resulting in a much thinner ^^F peak, so that the surface ^^F/^^0 
ratio is neither altered after the first, nor after the second dredge- up. 

Clearly, the first and second dredge-ups in intermediate-mass stars cannot account for 
even the lowest F overabundances observed in M giants. This discrepancy may indicate that: 

(i) the stellar fluorine abundances are systematically overestimated; or 

(ii) the solar-system fluorine abundance is in error; or 

(iii) the actual ^^F(p,a)^^0 reaction rate is slower than the intermediate rate provided by 
Kious (1990) and adopted in our calculations. In fact, that reaction rate is still largely 
uncertain in the temperature range of interest for hydrogen burning (see Arnould et al. 
1995). To investigate the impact of these uncertainties on the first dredge-up predictions, a 
new evolutionary sequence has been computed up to the first dredge- up for the 3Mq star, 
adopting the lower limit for the ^^F(p,Q!)^^0 rate. The surface fluorine overabundance then 
reaches [F/0] = 0.13 dex (see Fig. [l^ ). Although still on the low side of the observations, 
these predictions at least fall within the estimated ±0.1 dex uncertainty on the observational 
data. The case of a solar-metallicity low-mass (1.5 Mq) star has also been investigated with 
the low rate. The resulting surface F enhancement, also shown in Fig. |l^, is much smaller 
than in the 3M.Q star. 

Clearly, the elucidation of the differences between the solar-system fluorine abundance 
and that of normal K and M giants remains a challenge for future studies. In Sect. 5.2, 
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Figure 11: Fluorine abundance (thick solid line, corresponding to the left scale) in the 
envelope of the 6Mq, Z=0.02 star, and temperature (dotted line, to be read on the right scale, 
in units of 10^ K) at the base of the envelope. Time t = corresponds to the occurrence of 
the first convective pulse. The temperature drop after each pulse results from the expansion, 
and concomitant cooling, of the hydrogen-rich layers following the thermal pulses 




devoted to the analysis of the surface fluorine abundance resulting from the third dredge-up, 
we will therefore normalize our predictions in such a way that the pre-AGB ^^F abundance 
corresponds to the average value observed in solar- metallicity K and M giants (Fig. 14), 
whatever its origin might be. 



5.2 Third dredge-up 

After a pulse, the convective envelope extends inward and, in some cases, can penetrate the 
layers formerly involved in the thermal pulse. When this happens, the He-burning ashes, like 
^^C and ^^F, are mixed in the convective envelope and brought to the surface. This process, 
called third dredge-up, accounts for the gradual transformation of a M star into a S and then 
C star along the AGB. 

Though the very existence of S and C stars on the AGB supports the occurrence of such 
dredge-ups, our stellar models like most others fail in obtaining this process in a self-consistent 
way. No fully self-consistent comparison of the observed ^^F abundances with our predictions 
is thus possible. 

Nevertheless, a comparison is attempted by parametrizing the dredge-ups. Two parame- 
ters need to be specified for that purpose: (i) the depth A = AMd/{Ma + Mf,), where AM^ 
is the mass of matter from the pulse mixed into the convective envelope, and {Ma + Mf,) 
is the mass involved in the pulse (see Fig. ^); (ii) the pulse number ngtart after which the 
third dredge-up first occurs. It is then assumed that dredge-ups with the same depth A occur 
after each subsequent pulses. Though these free parameters introduce some uncertainties in 
our predictions, it was shown in Paper I (Eq. 1, Sect. 3.2) that the use of the (C/0, F/0) 
diagram strongly limits their impact. Actually, predictions almost independent of A can be 
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Figure 12: Mass fraction profiles of ^^C, ^^C, ^^N, ^^o, i^q, isq and ^^F before the first 
dredge-up in the M3Z02 star. The maximum extension of the convective core during central 
hydrogen burning (labeled core), and the maximum inward penetration of the envelope during 
the red giant phase (labeled envelope) are indicated at the bottom of the figure 



made by using this diagram. 

The parametric dredge-up calculations start with the envelope abundances at pulse n-start 
as given by the detailed models. These envelope abundances are then modified after each 
subsequent pulse by mixing into the envelope the matter from the underlying layers, until 
a fraction A of the ashes of the pulse has been dredged up. The abundance profiles are 
taken from the detailed model calculations. When HBB is taking place, as in the M6Z02 star 
(Sect. [4. 3D , its effect on the envelope abundances is incorporated in the parametric calculations 
by correcting the abundances by an amount equal to the relative variation obtained during 
the interpulse in the detailed calculations. 

Figure ^ displays how the surface (C/0, F/0) ratios vary in the M3Z02 and M6Z02 stars 
as a result of the dredge-ups. As discussed in Sect. 5T, our predictions are normalized in 
such a way that the surface F abundance prior to any dredge-up corresponds to the average F 
abundance observed in K and M giants. Predictions are given for ngtart = 1, 4 and 9, with a 
dredge-up depth A = 0.3. In the case of the M3Z02 star, the surface fluorine overabundances 
at a given ratio are the largest for ristart = 9 as expected, since the intershell fluorine 

abundance saturates from pulse 9 onwards in that star (Fig. However, even in this case, 
our predictions can only account for the lowest observed fluorine overabundances. Predictions 
with A = 0.1 lead to slopes in the (C/0, F/0) diagram similar to the ones obtained with 
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Figure 13: The diagram, comparing the abundances observed in several 

classes of red giants from Paper I with our predictions for the first dredge-up (small filled circle 
along the solid lines) and the second dredge- up (small filled square along the solid lines) . The 
thick lines correspond to the M3Z02 and M6Z02 stars, as labeled. The thin lines correspond 
to the first dredge-up predictions in 1.5 and SMq stars (italic labels) when the lower limit 
is taken for the ^^F(p,a)^^0 rate (Kious 1990), leading to the maximum ^^F production. 
Note that, unlike in Paper I, /^^O] = \og{^'^Y /^^0)^-\og{^'^Y /^^0)q, where subscripts * 
and refer to stellar and solar abundances, respectively. The error bar corresponds to the 
estimated uncertainty on the observational data 

A = 0.3, and are not reported in Fig. ^ for clarity. The above conclusion is thus independent 
of the dredge- up parameter A. 

In the M6Z02 star, the dredge-ups do not lead to any significant F overabundances. This 
is due to the combined effects of the large dilution of the intershell ^^F in the massive envelope, 
and to HBB. Even though a slight increase in the surface F abundance is predicted during 
the early pulses, due mainly to fluorine production in the envelope itself (see Sect. |4.3D , the 
operation of HBB at later pulses prevents any further increase. Moreover, by converting 
the envelope ^^C into ^^N, HBB leads to a decrease of the ratio in these massive 

AGB stars. Therefore, it is rather unlikely that the giants exhibiting large F/0 ratios have 
experienced HBB, unless their F/0 and C/0 ratios prior to the onset of HBB were even 
larger than their present values. But according to the present models, massive AGB stars 
are not very efficient in increasing their surface F abundance even in the absence of HBB. 

In this respect, WZ Cas is a puzzling case. This star is known to be super Li-rich (Denn, 
Luck & Lambert 1991, Abia et al. 1993), to have a low ^^C/^^C ratio, and to be enriched 
in s-process elements (Dominy 1985). Lithium is produced in the envelope of massive AGB 
stars (Sackmann & Boothroyd 1992) having temperatures Tg ^ 50 at their base (Mowlavi 
1995b). As no large surface F overabundances are expected in these typical HBB conditions, 
the combination of large Li and F excesses at the surface of WZ Cas remains unexplained so 
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Same as Fig. |l^, but for parametrized third dredge-up predictions. 



Figure 14: Same as Fig. |1^, but for parametrized third dredge-up predictions. Dotted 
hues correspond to predictions with A = 0.3 and n-start = 1, 4 or 9, as defined in the text. 
Predictions are normahzed in such a way that the surface F abundances prior to any third 
dredge-up correspond to the average F abundance observed in K and M giants. The soUd 
hne corresponds to the predictions for the M6Z02 star in the absence of third dredge-ups, 
reflecting the operation of HBB. Note that time is increasing from left to right (i.e. increasing 
C/0) in the M3Z02 tracks, but from right to left (i.e. decreasing C/0) in the M6Z02 ones. 
Underlined symbols denote stars with large N overabundances (see Paper I) 



far. More observations of F abundances in super Li-rich stars are needed to know whether 
such a combination is common among these stars. 

It is noteworthy that J stars (represented by open triangles in Fig. U) exhibit rather 
large F overabundances. This class of giants is characterized by low ^^C/^^C ratios and no 
s-process overabundances (Lambert et al. 1986, Dominy 1985). Since their N abundance is 
not especially high, HBB alone seems inadequate to account for their peculiar abundances. 
The combination of low ^^C/^^C and high ^^F/^^0 ratios is also puzzling in that respect, as 
the extensive hydrogen burning implied by the low ^■^C/^'^C ratio is also expected to destroy 
any pre-existing ^^F, unless hydrogen burning occurs at Tg < 20 (Sect. 4.3). 

In conclusion, the surface F abundances predicted for the intermediate-mass stars consid- 



ered in this paper are in rather poor agreement with the observations displayed in Fig. |14. 
The predicted F abundances in the M3Z02 star are just able to account for the lowest ob- 
served overabundances (at a given C/0 ratio). That conclusion is almost independent of 
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the parameters A and ngtart- Actually, a better match to the observations would require to 
increase the ^^F/^^C ratio in the pulses by a factor of at least 10 (see Eq. 1 and Fig. 9 of 
Paper I), i.e. X(^^F) ^ 210~^ in the pulse. Equation (^) indicates that the maximum ^^F 
production level in the pulse is fixed by the available ^^C supply. For secondary ^^C sup- 
plied by the CNO cycles operating in the hydrogen-burning shell, Xa(^^C) ~ 10~^, and the 
required ^^F production level is thus not achieved. That conclusion holds true for low-mass 
stars. It seems thus unavoidable that the most F-rich stars of Fig. 14 require a primary ^^C 
supply resulting possibly from a non-standard mixing of protons with the ^^C pocket left over 
by the pulse (e.g. Hollowell & Iben 1988). Since a similar need for a primary ^^C supply is 
in fact coming from the heavy-element overabundances observed in the same stars (Smith & 
Lambert 1990), one may expect the F and s-process overabundances to be correlated. The 
fact that the most F-rich stars are also the most heavy element-rich (Fig. 12 of Paper I) and 
the most N-rich (at a given C/0 ratio; underlined symbols in Fig. [l^ ) provides evidence for 
primary ^^C and ^''N supplies in these stars. 
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Figure 15: Intershell ^^O mass fractions after each pulse as a function of the maximum pulse 
temperature for the M3Z02 and M6Z02 stars. The dashed line corresponds to the surface 
^^O abundance of the M3Z02 star prior to any third dredge-up 

Another important constraint on the models comes from the surface ^^O abundance. In 
Paper II, large ^^O overabundances, not supported by the observations (Smith & Lambert 
1990), were predicted to result from the third dredge-up. The new ^'^O (a , 7) ^^Ne reaction 
rate (Giesen et al. 1994), which is faster than the CF88 rate used in Paper II by a factor of 
100 at Tg = 160, removes the discrepancy. Figure 15 displays the after-pulse ^^O intershell 
mass fraction as a function of the maximum pulse temperature, as well as the envelope 
abundance prior to any third dredge-up in the M3Z02 star. It is seen that the after-pulse 
^'^O abundance drops below the envelope abundance after the 4th pulse in the M3Z02 star. 
Therefore, in order to prevent surface ^^O overabundances, third dredge-ups should occur 
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Figure 16: Same as Fig. 14, but for M3Z001 (see text). The scales are the same as in Fig. ^ 



after the 4th pulse in that star. The ^^O overabundance may however be more difficult to 
avoid in low-mass stars, where pulses are cooler. 



5.3 Surface fluorine predictions in low-metallicity intermediate-mass stars 

Surface fluorine abundances resulting from parametrized third dredge-ups are predicted for 
the M3Z001 star in a similar way as for the solar-met allicity stars (Sect. |5.2| ). It is important to 
note that these predictions for the low-metallicity star rely on several important assumptions 
regarding the initial envelope abundances. In our models, all initial abundances are taken 
equal to the solar-system values, scaled to the metallicity by the factor Z/Zq. Such a scaling 
rule is justified for C and supposedly F, since [C/Fe] = in unevolved disk stars (Edvardsson 
et al. 1993) and [F/Fe] ~ is suggested by the observations of the low-metallicity K giant 
a Boo (Arcturus, [Fe/II]=-0.6; see Paper I). It however does not hold for O, since it is 
well known that [0/Fe] = —0.4 [Fe/H] in the range —1 < [Fe/H] < (e.g. Edvardsson et al. 
1993). The initial ^^O abundance in the envelope of the M3Z001 star has been corrected 
correspondingly. With these assumptions, the starting point of the dredge-up tracks in the 
(i2Q/i6Q^[i9p/i6o]) diagram corresponds to (0.12,-0.52) (see Fig. |l^). 

As can be seen on Fig. |l^, third dredge-up episodes in the M3Z001 star do not alter 
the surface fluorine abundance as strongly as in the solar-metallicity stars (Fig. U). The 
increase in [^^F/^^0] with respect to its initial value is only 0.1 dex after a dozen of pulses, 
to be compared with as much as 0.4 dex in the solar-metallicity stars. The less steep slopes 
predicted in the (C/0, F/0) diagram for the M3Z001 star result from the combination of two 
effects: (i) the lower fluorine production efficiency in the hot pulses of the low-metallicity star 
as compared to the ones in the solar-metallicity cases (compare Figs. ^ and P). The largest 
surface F overabundances in the M3Z001 star are therefore obtained when the dredge-ups 
occur after the early, cool pulses (compare tracks labeled 1 and 9 in Fig. 16); (ii) the primary 



nature of ^^C produced in the pulse, as compared to the secondary nature of ^^F in our 
models. This leads to ^^F/^^C ratios in the dredge-up material lower in low-metallicity stars 
than in solar-metallicity stars. 

In conclusion, rather small F overabundances are expected in low-metallicity AGB stars, 
in the case that the fluorine production results from secondary ^^C only. Observations of F 
in low-metallicity giants would deflnitely be of great interest. 
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6 Conclusions 



Detailed calculations of fluorine production have been performed in intermediate-mass stars 
with (M/Mq, Z) = (6, 0.02), (3, 0.02) and (3, 0.001). It is confirmed that ^^F is produced 
in pulses through ^^N (a , 7) ^^F (P~^) ^^O (p , a) ^^N (a , 7) -^^F , the necessary protons coming 
from (n,p) reactions with neutrons supplied by ^^C (a , n) ^^O . Three different regimes of 
fluorine production can be identified: 

1) In pulses with maximum base temperatures Tg ^ 260, ^^N is produced within the pulse. 
For Tq ^ 220, ^^N (a , 7) ^^F is however too slow to produce significant amounts of F. The 
maximum fluorine production is reached for 220 ^ Tq ^ 260, at a level limited by the ^^C 
supply and the efficiency g of ^^N production in the reaction chain starting at -"^^C. Typical 
values for the efficiency are g ~ 0.2, and depend mainly upon the fraction of available neutrons 
involved in (n,p) reactions. ^^N and ^^Al are the main neutron-to-proton converters. The 
efficiency is found to depend also on the overlap factor between successive pulses. 

2) In pulses with Tq ^ 260, the efficiency of -"^^N production within the pulse drops dramat- 
ically. However, the temperature of the intershell region now gets high enough for ^^N to be 
produced in radiative layers during the interpulse period. When injected in the next pulse, 
that ^^N is entirely converted into ^^F by a-captures. The efficiency of ^^N production is lower 
in radiative layers than in a pulse {grad ~ 0.1). The intershell ^^C nuclei not burned during 
the interpulse period are ingested by the pulse close to the time of maximum pulse temper- 
ature. At these high temperatures, ^^N(a,7) "^^F is fast enough to substantially reduce the 
number of ^^N seeds available to ^'^N (n , p) ^^C . The overall efficiency of ^^F production is 
thus lower in this regime than in regime 1. 

3) In pulses with Tq ^ 300, ^^F (a , p) ^^l^e destroys ^^F. 

The evolution of the after-pulse fluorine abundance Xout(^^F) in the intershell zone is 
dictated by the succession of these three regimes. It first increases during the first regime until 
a maximum value is reached, given by XautC^F) = §gX{^^C) (where X{^^C) corresponds 
to the mass fraction of the available ^^C supply). During regime 2, XontC^F) drops until 
it reaches a plateau, with a level fixed by the efficiency grad- Finally, Xout(^^F) decreases 
steadily when regime 3 operates in the late pulses. The optimum F production in AGB stars is 
achieved for thermal pulses having maximum base temperatures in the range 220 ^ Tq ^ 260. 

All three regimes have been encountered in the M3Z001 star where the pulse temperatures 
reach Tq = 260 at the 4th pulse and 300 at the 13th. The plateau corresponding to the second 
regime is reached after 9 pulses in that star. The M6Z02 star enters regime 2 at the 7th pulse, 
while the M3Z02 star is still in regime 1 after 13 pulses. 

Comparing our predictions of surface fluorine abundances with the observations leads to 
the following conclusions: 

i) Neither the first nor the second dredge-ups can account for the observed discrepancy 
between the solar-system fluorine abundance and that of normal K and M giants. 

ii) Our predictions of surface fluorine abundances resulting from parametrized third dredge- 
ups in intermediate-mass stars can only account for the lowest F overabundances observed 
in AGB stars. That conclusion is expected to hold true for low-mass stars when fluorine pro- 
duction results from secondary ^^C only. No concomitant ^^O overabundances are predicted 
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in our intcrmcdiatc-mass stars, in agreement with the observations. 

iii) The largest F overabundances observed in C stars call for an additional source of ^^C, 
of primary origin. A larger ^^C supply could account for the observed overabundances of 
both F and s-process elements. Detailed computations are however necessary to confirm that 
prediction. 

iv) According to our models, massive AGB stars are not expected to build up large surface 
F abundances. Moreover, when operating at Tq > 70, HBB leads to a gradual destruction of 
the envelope fluorine, thereby limiting the surface F overabundance that would result from 
the third dredge-up. The large fluorine overabundance reported for the super Li-rich star 
WZ Gas (where HBB is supposed to be operating) remains therefore unexplained so far. 

v) Hot thermal pulses in intermediate-mass, low-metallicity AGB stars are less efficient in 
producing F. The lower ^^F/^^C ratio from these pulses does not lead to large F overabun- 
dances at the surface of these AGB stars. Any large overabundances that would be reported 
by (future) observations of low-metallicity stars, would require a primary source of ^^G. 
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Appendix: references to nuclear reaction rates 

This appendix provides the references for the reaction rates used in the present study when 
different from CF88 or from Beer et al. (1992, for neutron-capture reactions). Tables Al, A2 
and A3 refer to proton-, neutron- and ct-capture reactions, respectively. The Hauser-Feshbach 
method, implemented in the SMOKER code (Thielemann, Arnould &: Truran 1987), has been 
used to derive the reaction rates for proton captures on ^^Si, ^^Si, ^^P, ^^P, ^^P, ^^S, ^^S and 
^"^S, not listed in Table Al. 
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Table Al: References for proton-capture reaction rates, if different from CF88 



lie (p,7) 12n(/3+)12c 


Descouvemont & Baraffe 1990 


(p,7) (/3+)14n 


Decrock et al. 1991 


i^C (p , 7) i-'^N 


Wicscher et al. 1990 


1^0 (p , a) i^N 


Berheide et al. 1993" 


i70(p,7)i«F 


Berheide et al. 1993" 


i^F (p , a) 


Kious 1990 


2iNe(p,7) ^^Na 


Gorres et al. 1982, Gorres 




et al. 1983 


22Ne (p,7) 23Na 


Gorres et al. 1983 


22Na(p,7) ^^Mg (/3+)23Na 


Scuthc ct al. 1990^ 


23Na (p , 7) 24Mg 


CF88 + G6rres et al. 1989^= 


25Mg (p,7) 26Alg 


Iliadis et al. 1990^^ 


25Mg (p,7) 


Iliadis et al. 1990*^ 


26Mg (p,7) 27A1 


Iliadis et al. 1990^ 


(p,7) 27Si(/3+)27Al 


Champagne et al. 1993" 


27Al(p,7) ''Si 


Timmermann et al. 1988" 


(p , a) 24Mg 


Timmermann et al. 1988 




+ Champagne et al. 1988^'" 


28Si(p,7) (/?+)29Si 


Graff ct al. 1990 


(p,7)33Cl(/?+)33S 


Ihadis et al. 1992 



" Geometric average between lower and upper limits 

Resonance at 70 keV from Berheide et al. (1993), plugged in expression from Landre et al. 
(1989) with fi=f2=l 

Resonant contributions from Gorres ct al. (1989) 

Geometric average between lower and upper limits for Tg > 0.01 
^ Resonant contributions from Champagne et al. (1988) 
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Table A2: References for neutron-capture reaction rates, if different from Beer, Voss & 
Winters (1992) 



i^C (n , 7) 




Ohsaki et al. 1994 


"C (n , 7) 




Raman et al. 1990 


i^N (n , p) "C 




CF88 (reverse) 


i^N (n , p) 




Koehler h O'Brien 1989 


(n , 7) °N [p ^ 


IDQ 


FCZ67 


^^0 (n , p) ^^N 




CF88 (reverse) 


^'^0 (11,7) ^^0 




Wagoner 1969 


I'O (n,a) 




Schatz et al. 1993 


^*^0 (n , 7) ^^0 [fi 


1 Q-rn 

) ^^F 


Rauscher et al. 1994 


^'^F (n , a) ^''N 




Thibaud 1989 


(n , p) 




Thibaud 1989 


2iNe (n , 7) 22Ne 




Almeida & Kappeler 1983 


26A1S (n,7) 27A1 




WFHZ78 


26AF(n,p) 26Mg 




CF88 


26Alg (n , a) 23Na 




CF88 


26AP (n,7)27Al 




WFHZ78 


26AP (n,a) 23Na 




CF88 


26AP (n,p) 26Mg 




CF88 


3iSi(n,7)32Si 




SMOKER 


32Si (n,7) 33Si (/?- 


)33p 


SMOKER 


32p (n,7)33p 




SMOKER 


32p(n,p)32Si 




SMOKER 


33p (n,/3-)34s 




SMOKER 


33s (n,p)33p 




SMOKER 


33s (n,a) 30Si 




Wagemans et al. 1987 



Table A3: References for a-capture reaction rates, if different from CF88 



[a 


,7)^^0 


Caughlan ct al. 1985 


13c (a 


, n) 160 


Drotleff et al. 1993 


i^C (a 




Hashimoto et al. 1986, Gai et 






al. 1987, Funck h Langanke 






1989 


i^N (a 


,7)i9F 


de Oliveira et al. 1995 


1^0 (a 


,7)22Ne 


Giesen et al. 1994 


22Ne (a , n) 25Mg 


Drotleff et al. 1993 


22Ne (a, 7) 26Mg 


Drotleff et al. 1992, priv. 
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